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Ion H2 can dissociate in a strong magnetic 
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Abstract 

In framework of a variational method the molecular ion H% in 
a magnetic field is studied. An optimal form of the vector poten- 
tial corresponding to a given magnetic field (gauge fixing) is cho- 
sen variationally. It is shown that for any magnetic field strength 
as well as for any orientation of the molecular axis the system 
(ppe) possesses a minimum in the potential energy. The stable 
configuration always corresponds to elongation along the mag- 
netic line. However, for magnetic fields B > 5 x 10 11 G and some 
orientations the ion becomes unstable decaying to if-atom+p. 



Ion H2 is the simplest one-electron molecular system, which is more sta- 
ble than hydrogen atom. It appears to be one of the most studied problems 
in non-relativistic quantum mechanics. In particular, the wealth of physical 
phenomena displayed by this system when placed into a magnetic field be- 
comes of a great importance in astrophysics, solid state and plasma physics. 
For instance, as the magnetic field grows the system becomes more and more 
strongly bound and compact. Such a behavior leads naturally to a guess 
that in spite of the huge temperatures on neutron star surfaces their at- 
mosphere can still contain molecular objects |l], [|. On the other hand, a 
shrinking of the equilibrium distance between protons with the growth of 
magnetic field, increases drastically the probability of nuclear fusion |J. It 
is quite surprising that such a shrinking is also accompanied by a change 
from ionic to covalent character at~5-10 11 G[|J (see also ||). The goal of 
this Letter is two-fold. Firstly, to show that the system (ppe) has always a 
minimum and correspondingly the molecular ion can exist for magnetic 
fields < 4.414 ■ 10 13 G. Secondly, to demonstrate that for B > 5 x 10 11 G 
and for some orientations of the molecular axis the ion becomes unstable 
dissociating to if-atom + p. The variational method is used to study this 
problem. 

The Hamiltonian which describes the molecular ion placed in a uni- 
form constant magnetic field directed along the z— axis, B = (0, 0, B) (see, 
for example, 0) is given by 

H = p* + l---- + (pA) + A\ (1) 
R r\ r 2 

(see Fig.l), where p = — zV is the momentum. A is a vector potential, which 
corresponds to the magnetic field B. 

The vector potential for given magnetic field is defined ambiguously, up 
to a gauge factor. Thus, the Hamiltonian is gauge-dependent but not the 
observables. Since we are going to use an approximate method for solving ([!]) 
our energies can be gauge- dependent (only the exact ones would be gauge- 
independent). Hence one can choose the form of the vector potential in 
an optimal way. Let us consider a certain one-parameter family of vector 
potentials corresponding to the constant magnetic field B (see for example 

i) 

A=B(-(l-£)y,Zx,0) , (2) 

where £ is the parameter to be fixed in a certain optimal way. If £ = 1/2 
we get a gauge called symmetric or circular, while £ = corresponds to an 
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asymmetric gauge (see |J). By substituting (^) into (|T]) we arrive at the 
Hamiltonian 

H = -V 2 + ^----+iB[-(l-0yd x +Cxd y }+B 2 [ex 2 + (l-0V] ■ (3) 
K ri r 2 

The idea of choosing an optimal gauge has widely been exploited in quan- 
tum field theoretical considerations. It has also been discussed in connection 
with the problem at hand (see for instance and references therein). Per- 
haps, the first constructive (and remarkable) attempt to apply this idea was 
made by Larsen In his study of the ground state it was explicitly shown 
that gauge dependence of energy can be quite significant and even an over- 
simplified optimization procedure improves the numerical results. 

Our aim is to study the ground state of (|3]). It is not difficult to see 
that there exists a certain gauge for which Hamiltonian (|3|) has a real ground 
state eigenfunction. This gauge will be here sought after and correspondingly 
deal with real trial functions in our variational calculations. In this case 
one can prove that the expectation value of the term ~ B in (|3p vanishes 
when it is taken over any real normalizable function. So, without loss of 
generality we can omit this term in the Hamiltonian. Finally, the recipe of 
our variational study can be formulated as follows: Construct an adequate 
variational real trial function J77]/ ; which reproduces the original potential 
near Coulomb singularities and at large distances, where £ should be included 
as a parameter. Perform a minimization of the energy functional by treating 
the trial function's free parameters and £ on the same footing. In particular, 
such an approach enables one to eventually find the optimal form of the 
Hamiltonian. The above recipe was successfully applied in a previous study 
of in magnetic field || and led to predict the existence of the exotic ion 
/tf+atB^lO 11 *?!. 

One of the simplest trial functions satisfying the above-mentioned crite- 
rion is 

_ e -ai(ri+r 2 )-B[/3i^ 2 +Ay(l-0?/ 2 ] ? (4) 

(cf. 0) where a%, Pi x ,iy are variational parameters. We here assume that 
£ G [0, 1] which is a restriction that will later be justified. Actually, this is a 
Heitler-London function multiplied by the lowest Landau orbital associated 
with the gauge (Q). Presumably this function describes internuclear distances 
near the equilibrium and a covalent character. Another possible trial function 
is 

v|> 2 — ^e~ a2ri + e~ a2T2 ^j e - B ^& 2 +P2y{i-£,)y 2 } ^ ^ 
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(cf. H) where a 2 , (3 2x ,2y are variational parameters. This is a Hund-Mulliken 
function multiplied by the lowest Landau orbital. One can assume that for a 
sufficiently large internuclear distance R this function dominates, thus giving 
an essential contribution in this regime. Hence, it describes an interaction of 
a hydrogen atom and a proton (charged center), and can also describe a pos- 
sible decay mode of H£ onto them. There are two natural ways— linear and 
non-linear— to incorporate the behavior of the system both near equilibrium 
and at long distances in a single trial function. The non-linear interpolation 
is of the form 

$ 3 _ 1 = f e - a 3ri-a4r2 _|_ e -a 3 r 2 -a i r 1 e -B[f3 [ix £x 2 +f3 [iy (l-£)y 2 ] ^ (g) 

(cf. ||) where a^^, (3^ x ^ y are variational parameters. This is a Guillemin- 
Zener function multiplied by the lowest Landau orbital. The linear superpo- 
sition is given by 

* 3 -2 = A x ¥ x + A 2 V 2 , (7) 

where one of the parameters A X}2 is kept fixed. The final form of the trial 
function is a linear superposition of functions @ and (0) 

% rial = + A 2 ^ 2 + A 3 _ 1 ^ 3 _ 1 , (8) 

where only two out of the three parameters A's are variationally treated. The 
total number of variational parameters in (g) is fourteen when £ is included. 

It is easy to prove a general statement that if a system possesses axial 
rotational symmetry (in our case it appears if the molecular axis coincides 
with the magnetic line, 9 = 0°, see Fig. 1), the optimal gauge corresponds 
to £ = 1/2 (symmetric or circular gauge). It is precisely the gauge which 
was used in most of previously performed H 2 -studies. However, this is not 
the case if 9 7^ 0°. As an example one can see in Fig. 2 the behavior of 
£ as a function of 9 at B = 10 12 G. It is typical behavior for all studied 
magnetic fields both weak and strong, up to the non-relativistic limit, B = 
4.414 x 10 13 G. It justifies our above assumption with regard to the domain 
of gauge parameter, £ e [0, 1]. 

We carried out extensive studies of the ground state la g of H 2 for mag- 
netic fields B = — 4.414 x 10 13 G and orientations ranging from 0° (parallel 
configuration) to 90° (perpendicular configuration). They turned out to be 
more accurate than any available results yet obtained except for a domain of 
small magnetic fields for the perpendicular configuration, where the results 
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by Wille |TT| appear to be slightly better Q The detailed numerical analysis 
and comparison with available calculations will be published elsewhere. 

As previously obtained by other authors |3|, |TTj, [|] we confirm quantita- 
tively a natural expectation that the parallel configuration is the most stable 
for all magnetic fields B < 4.414 x 10 13 G, where non-relativisitic consider- 
ations are valid. The total energy of the molecular ion H2 as a function 
of angle 6 for different magnetic fields, is shown in Fig. 3. For any mag- 
netic field in the region B = — 4.414 x 10 13 G and any orientation a well- 
pronounced minimum in the total energy is attained at finite internuclear 
distances. This is in contradiction with a statement by Khersonsky H about 
the non-existence of a minimum for some values of the magnetic field and ori- 
entation. Perhaps, it it worth to emphasize that in that article the variational 
study was carried out using a trial function which is almost coincident with 
that of eq.(^). We can thus presume that the above statement is an artifact 
arising from insufficient accuracy of calculations. The horizontal line in Fig. 
3 presents the hydrogen atom total energy in the magnetic field (see ||). For 
magnetic fields B > 1.8 x 10 11 G the total energy of the atom becomes lower 
than that of H£ for angles larger than some critical angle, 8 cr . It points to 
the possible dissociation channel — > H-atom + p. The dependence of the 
critical angle 8 cr on the magnetic field is shown in Fig. 4. It is quite striking 
that dissociation occurs for a wider and wider range of orientations as the 
magnetic field grows reaching 25° < 6 ^ 90° for B = 4.414 x 10 13 G. 
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1 These results were obtained using a basis set expansion including up to approximately 
500 terms depending on the value of magnetic field. 
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Figure 2: Gauge parameter as a function of the inclination angle for H^- 
The magnetic field B = 10 12 G was taken as an example. This dependence is 
found in the present study. 
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Figure 3: Total energy of the ion as a function of the inclination ang 
The horizontal lines refer to the H ground state energy taken from Lai et 
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